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Abstract 



We study the three-body system with short-range interactions charac- 
terized by an unnaturally large two-body scattering length. We show that 
the off-shell scattering amplitude is cutoff independent up to power correc- 
tions. This allows us to derive an exact renormalization group equation for 
the three-body force. We also obtain a renormalized equation for the off-shell 
scattering amplitude. This equation is invariant under discrete scale transfor- 
mations. The periodicity of the spectrum of bound states originally observed 
by Efimov is a consequence of this symmetry. The functional dependence 
of the three-body scattering length on the two-body scattering length can 
be obtained analytically using the asymptotic solution to the integral equa- 
tion. An analogous formula for the three-body recombination coefficient is 
also obtained. 
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Recently, there has been renewed interest in the three-body system with short-range 
interactions [1-9]. From the effective field theory (EFT) perspective, an understanding of 
this system is an important ingredient for a successful description of few- and many-body 
systems in nuclear and atomic physics [10-13]. Despite its simplicity, this system exhibits 
many interesting features [14-16]. For example, when the two-body scattering length a2 
approaches infinity, the three-body system exhibits an infinite number of shallow bound 
states [15]. The equations describing three particles interacting via strong short-range two- 
body forces have been known for a long time [17]. While these equations are well behaved 
in some channels, they do not have a unique solution for spinless bosons or three nucleons 
with total spin ,5 = 1/2 [18,19]. 

The renormalization of the EFT for the three-body system with short-range interactions 
has been discussed in detail in Ref. [3]. Here we give a brief summary of the results. 
For nonrelativistic particles interacting via short-range forces, the Lagrangian consists of a 
nonrelativistic kinetic term and an infinite number of contact interactions with an increasing 
number of derivatives. For systems in which the two-body scattering length 02 is much larger 
than the characteristic range of the interactions, the leading two-body contact interaction 
with no derivatives needs to be treated nonperturbatively [20-22]. For three-body systems, 
the leading order Feynman diagrams can be summed using an integral equation which is 
identical to that of Ref. [17] if the three-body force is of natural size and subleading. For 
three bosons or for the spin- 1/2 state of three nucleons, this integral equation exhibits strong 
dependence on the cutoff used to regulate the theory. From the viewpoint of effective field 
theory, this is surprising as the Feynman diagrams that the integral equation is designed to 
sum are individually finite. Their sum, however, does not converge and this is the origin of 
the cutoff dependence in the integral equation. 

The cutoff dependence of the three-body integral equation is properly interpreted via 
renormalization theory. In this case, the sensitivity to the cutoff indicates that a three-body 
contact interaction, which one would regard as subleading on the basis of naive dimensional 
analysis, is in fact leading order. In Ref. [3], it was shown that the cutoff dependence 
in the three-body equation can be properly renormalized with the inclusion of this three- 
body force. The three-body force exhibits a very unusual renormalization group flow: it 
is characterized by a limit cycle. The relevance of a single three-body operator provides a 
compelling explanation for the existence of the Phillips line [16]. The effective theory has 
enjoyed successful phenomenological applications in neutron-deuteron scattering [1,2,4-6], 
the scattering of ^He molecules [3], and three-body recombination of atoms in Bose- Einstein 
condensates [7]. 

In this paper, we will show that once the three-body force is included, the cutoff depen- 
dence of off-shell three-body amplitudes vanishes as the cutoff is taken to infinity. This fact 
allows us to write down renormalized equations in which the cutoff is completely removed. 
These renormalized equations exhibit invariance under discrete scale transformations. This 
invariance is exact at leading order in the effective theory and can be used to constrain 
the functional form of three-body observables. The spacing of energy levels of low-lying 
three-body bound states is a direct consequence of this symmetry. The invariance of the 
three-body observables under the discrete scale transformations and the spacing of the en- 
ergy levels has been derived previously in a very different manner by Efimov [15,23]. We 
also show how asymptotic solutions to the equations can be used to determine analytically 



the dependence of the three-body scattering length and the three-body recombination rate 
on the two-body scattering length. The formula for the three-body scattering length has 
also been derived in earlier work by Efimov [15,23], while the behavior of the three-body 
recombination rate has been extracted from fits to numerical solutions of the equations [7] . 

We begin with the integral equation for the elastic scattering of a particle and a bound 
state of the other two particles: 

[F(p; k)]-^K{k, p) = M{k, p; k;A) + ^ £ dq M{q, p; k; A)p(-^^) K{k, q) , (1) 

where 




F{p- k) 
M{q,p;k;A) 



Equation (1) has been derived in Ref. [3]. Here 02 is the two-body scattering length, and the 
total energy is mEi. = Sk^/A — l/a|. K{k,p) describes the S-wave scattering of a particle 
and a bound state with momentum ±k into a state with momentum ±p. This is an off-shell 
quantity except atp = k, where it is related to the scattering phase shift, k cot 6 = 1/K{k, k). 
The three-body scattering length is simply as = —K{0, 0). The integral equation is regulated 
by the ultraviolet cutoff A and H{K)/h? is the contribution from the three-body force. H{K) 
depends on the three-body parameter. A*, and evolves in such a way as to render the solution 
insensitive to A. 

We concentrate now on threshold scattering. Setting A; = and defining K{Q,p) = a{p), 
we obtain the equation: 



[Fip;0)]-'a{p) 



1 



p 



H(A) 2 

A2 TT 



dqM{q,p;0;A) a{q) . 



(2) 



It is clear from the numerical solutions of Eq. (2) in Ref. [3] that a{p) is independent of A up 
to power suppressed corrections. This fact is crucial for the derivation of the renormalized 
equation. On-shell quantities such K{k, k) are necessarily cutoff independent if the theory is 
properly renormalized. However, a{p) is an off-shell quantity for p 7^ 0, so it is not obvious 
that a{p) should be cutoff independent as well. To see why this is the case we will derive a 
renormalization group equation (RGE) for a{p). We begin by noting that it is possible to 
derive an equation in which the three-body force is eliminated by taking the difference of 
Eq. (2) with p ^ and the same equation with p = 0, 
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q"^ + qp + p"^ + 1/a: 



q'^ — qp + p'^ + l/al 



a{q) 



Acting on Eq. (3) with A d/dA, and using the fact that a(0) must be cutoff independent, we 
derive the following RGE for a{p): 
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a(A) 
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(4) 



Aj^a(g). 



The first term on the right hand side of Eq. (4) is of order a(A) Q^/A^, where A » Q = 
p, 1/02- In the hmit of large A, Eq. (4) becomes a homogeneous integral equation: 



[F{p;Q)]-'A-^a{p) = ^£dq 
which is trivially solved by 



2pq 



In 



q'^ + qp + p'^ + l/al 



q'^ - qp + p'^ + l/al 



g2 + l/al 



aA,(,) = o. 



(5) 



(6) 



Thus Eq. (5) provides an explanation for the cutoff independence of a{p). It is conceivable 
that Eq. (5) has other solutions, however, the numerical calculations of Ref. [3] (cf. Fig. 6) 
show that Eq. (6) is the correct solution up to corrections that vanish as A goes to infinity. 
We have not yet obtained an analytic proof of Eq. (6). 

It is possible to repeat this analysis for the off-shell amplitude K{k,p). After eliminating 
H{A.) from Eq. (1), acting with Ad/dA, and dropping power suppressed terms, we find 



[F{p;k)]-'A-^K{k,p) = [Fip;0)]-'A^K{k,0) 

mEu 



+ - / dq 

TT Jo 
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(7) 



Note that K{k,p) is not symmetric so K{k, 0) 7^ a{k). Eq. (7) is solved by 

A-^K{k,p)^Q. 



(8) 



In Fig. 1, we have plotted ^(0.8/02,]?) for four different values of the cutoff and a2A* = 16. 
The function is cutoff independent up to corrections of 0{l/{Aa2)). For the lowest value of 
the cutoff, Aa2 = 20, these corrections can be seen in Fig. 1. For the larger values of A, they 
are negligible. 

We can now use the cutoff independence of a{p) to derive the RGE for the three-body 
force. Setting p = in Eq. (2) gives the following equation for the three-body scattering 
length 03: 



3a2 2 HiK) 2 
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A2 



TT Jo 
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(9) 



where we have used a(0) = —03. All cutoff dependence in this expression is explicit. Taking 
derivatives with respect to A we obtain the exact RGE 
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FIG. 1. The function K{k,p) at ka2 = 0.8 for a2A* = 16 and four different values of the cutoff A. 



Prom a study of the asymptotics of Eq. (2) (cf. Refs. [18,3]), we know the behavior of a{p) 
for large p: 



a{p) = a2C{a2K) cos(so ln(p/A*)) , 



(11) 



where Sq ~ 1.006 and A* is the three-body force parameter. C is an unknown dimensionless 
function of 02 A* and expected to be of 0(1). We can substitute this asymptotic solution for 
a{p) into Eq. (10) in order to obtain an approximate RGE for H{A). If we neglect terms of 
C(l/(Aa2)) we find 



A 



_d_ ( HjAY 
dK [ A2 



-(1 + sg) cos(goln(A/A,)) 
cos(soln(A/A*) + sq sin(so ln(A/A*)) 




(12) 



which is solved by 



sin(soln(A/A*) — arctan(l/so)) 
sin(soln(A/A^,) + arctan(l/so)) 



(13) 



This solution for H{A) was previously obtained in Ref. [3] by requiring the low-energy 
solution a{p ~ 1/02) to be invariant under finite changes of the cutoff. Equations (12, 13) 
receive corrections that scale as l/(a2A), which could in principle be computed with the 
help of Eq. (10). In practice, Eq. (13) is an excellent approximation to the exact evolution 
ofH{A) (cf. Fig. 8 of Ref. [3]). 

We are now in a position to write renormalized equations for K{k,p) and a{p). Cutoff 
independence of K{k,p) implies that we are free to choose the cutoff to be whatever we like. 
Up to corrections of 0(l/(a2A)), H{A) vanishes if A = A„ where 
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A„ = A* exp 



1 / /I 
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(14) 



Setting A = A„ in Eqs. (1,2) results in 
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(16) 



These are interesting equations because all dependence on A has been removed in favor of 
the physical parameter characterizing the thrcc-body force, A^,, which appears in the upper 
limit of the integral. These equations receive corrections which fall off as powers of A. The 
leading corrections in 1/A come from the running of H[K) and scale as l/(a2A). They are 
suppressed in the limit a2 — > oo and can be made negligible by choosing n sufficiently large. 
Note that we are not taking the cutoff to infinity in a continuous manner, rather in a series 
of discrete steps for which H{A) = 0. 

The renormalizcd equations have an exact symmetry that is a consequence of the freedom 
to choose an arbitrary integer value for n in Eq. (14). This corresponds to rescaling A* 
by a factor of exp(n7r/so) and holding all other dimensionful quantities fixed. Note that 
this symmetry is trivially satisfied in the two-body system. As a consequence, physical 
observables should not change if the three-body force parameter is rescaled by exp{nn/ Sq). 

This symmetry can also be seen in the unrenormalized equation. The amplitude a{p) is 
a function of three variables: a{p] a2, A.^)- This function has a remnant of scale invariance 
due to the periodicity of H{A). If we perform the scale transformations: 



P 

l/a2 
A 

a(p;a2. A*) 



ap, 
a/a2 , 
aA , 

a~^a{ap; a~^a2, A* 



(17) 



in Eq. (2) where the cutoff is still present, we find 



[F(p; 0)] ^a{ap;a ^a2,A 
1 



(18) 



p 



^ + l/al 



H(aA) 2 
+ \ . ^ + 



A2 



2 Z'^'- 

— / dqM(q,p;0;aA)a(aq;a~^a2,A^). 
TT Jo 



Note that all dimensionful quantities except for A* are rescaled in the transformation in 
Eq. (17). Up to corrections of 0{l/{a2A)) (which can be made arbitrarily small by choosing 
an appropriate cutoff A), the three-body force runs according to Eq. (13) and has a hmit cycle 
with period exp{n7i / Sq) . As a consequence, if a(p; a2. A*) is a solution, then a{ap, a^^a2, A*) 
is a solution as well for all a — exp(n7r/so)- 



6 



This symmetry leads immediately to the unique spectrum of thrcc-body bound states 
originally discovered by Efimov [15]. The equation for the thrcc-body bound state can easily 
be obtained from Eq. (2). Since the bound state solutions correspond to standing waves, 
we drop the inhomogeneous term and replace with — ^3. It is sufficient to consider the 
equation in the limit k — 0, and we have 

where a{p) is related to the bound state wave function. The binding energies are then given 
by those values of for which Eq. (19) has nontrivial solutions. Equation (19) has the same 
symmetry as Eq. (2) if we transform — > a^B^. Consequently, if Eq. (19) has a solution 
for one value of B3 ~ l/(ma|), the symmetry then implies the existence of a spectrum of 
bound states with binding energies [15] 

B3" = i(..ln(aA))^^Pe^, (20) 

where n > is an integer. The prefactor d{so\n{a2A^)) is a periodic function of order one 
and has been calculated numerically in Ref. [3]. Because a stable three-body bound state 
cannot have a binding energy smaller than the two-body binding energy B2 — l/(ma|), n 
has to be larger than zero. Furthermore, due to the presence of the momentum cutoff A 
the maximal binding energy is of order A^/m. Inverting Eq. (20), we then find for the total 
number of bound states [15] 

N^-\n{Aa2). (21) 

TT 

If we increase the cutoff, the maximal binding energy increases as well [14]. The shallow 
bound states, however, are kept at constant binding energy through the running of H with 
the cutoff A [3]. Corrections to Eq. (20) from irrelevant operators have been studied in 
Ref. [9]. 

The discrete scale symmetry plays an important role in constraining the functional de- 
pendence of other three-body observables. We will consider the three-body scattering length 
03, the three-body effective range parameter r^, and the recombination coefficient a. In the 
recombination process, three atoms start at rest, then two of the atoms form a bound state, 
and the third atom balances energy and momentum. The recombination rate for three atoms 
in a gas of density n is u = an^, where a is the matrix element squared for the transition to 
take place. The integral equation that determines a is very similar to that of a(p) and will 
be given below. The quantities as, rs and a depend only on 02 and A* and the requirement 
that observables be invariant under A* — > exp(7r/so)A* highly restricts their functional form: 

a3(a2. A*) = a2 /(sq ln(a2A*)) , (22) 
r3(a2. A*) = 02 9^(50 ln(a2A*)) , 

Q;(a2,A*) = — /i(soln(a2A*)) , 

where f{x),g{x) and h{x) are all periodic functions: f{x) — f{x + 7r). The dependence on 
factors of 02 and m is determined by dimensional analysis. 
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In the following, we will investigate the observables in Eq. (22) in more detail. It turns 
out that the functional form of and a can be obtained from the asymptotic solution. We 
first turn to the three-body recombination rate a. For three atoms approximately at rest in 
the initial state, the final momentum for the bound state is p/ = 2/[\^a2). In Ref. [7], it 
was shown that a — where h{p) is the solution of 



327rV2 / 1 H{A) 
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TT io 
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dqq%{q) 
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It is straightforward to show that ph{p) satisfies the same asymptotic equation as a{p) from 
Eqs. (2, 16). Consequently, the ultraviolet behavior of a{p) and pb{p) is the same and 
Eq. (23) is renormalized by the same three-body force [7]. Furthermore, b{p) is independent 
of A up to power corrections, and we can choose the cutoff as in Eq. (14) to obtain a 
renormalized equation. 

The asymptotic solution for b{p) is expected to be valid only for p ^ 1/0,2. However, 
in the numerical calculations of Ref. [3] it was observed that the asymptotic solution works 
surprisingly well down to momenta of order l/a2. In this case we should be able to use the 
asymptotic solution to get the functional form of a. The solution to the asymptotic equation 
for b{p) has the form: 



bip) 



a2 cos(soln(p/A*)) 



P 



m 



(24) 



where the factors of 02 and m are determined from dimensional analysis. Note that the 
asymptotic equation does not determine the normalization of the solution. Evaluating the 
formula for a, we obtain 



« = \KPf)\^ 0^ — cos^[soln(a2A*) - 0.1] 



m 



(25) 



The recombination coefficient was computed numerically in Ref. [7]. The numerical solution 
can be fit with the following expression: 



a 



m 



68cos^[soln(a2A*) + 1.7] . 



(26) 



Equation (26) is invariant under A* exp(7r/so)A* as expected. The asymptotic solution 
gives the correct functional form but the normalization is not predicted. It does not predict 
the correct phase in Eq. (26) because the phase of the full solution at p ~ l/a2 is not equal 
to the asymptotic phase. 

Now we turn to the three-body scattering length which is given by 03 = — a(0). For 
p < l/a2, the asymptotic form is not a good approximation anymore. In fact, it does not 
even have a well defined limit as p — >^ 0. The only extra piece of information we need is the 
existence of momenta for which the value of a{p) is independent of A*. These meeting points 
were first observed in Ref. [3] (cf. Fig. 5). The first of these points occurs at po = l-l/a2, 
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and additional points occur at p„ = po cxp^nir / sq) with n > an integer. The existence and 
location of these points is not determined by the asymptotic equation, but from matching 
to the solution for p ~ 1/02- While the meeting points occur at off-shell momenta and are 
not observable, they severely constrain the form of the solution for p > 1/02- Taking into 
account the meeting points, the solution for momenta p ^ l/a2 is [3] 

/ X ^/ cos(soln(p/A,)) 

a{p) ~ 02 C — -- — T-—- , (27) 

cos(so ln(po/A*)) 

where C is a constant that is independent of 02^*. 

We now insert this form into Eq. (16), set p = 0, and keep only the terms that dominate 
in the limit 02 — > 00. We split the integral over q into an integral from to po and an 
integral from pg to A„. For large 02 the first integral is parametrically suppressed and can 
be neglected. The second integral can be performed by expanding the kernel in a formal 
power series, 

02 3 37ripo ~[ cos(soln(po/A*)) 

Ifi 39 °° ( — ^\^+^(n n U-2n 

= -y - 1;^^' E ,2^ (2^-1)^ [^^^ " + '° " • ^^^^ 

In Eq. (28), we have kept only the leading terms for for large 02. The three-body scattering 
length obtained from a fit to numerical solutions in Ref. [7] is 

aa « 02(1.4 - 1.8 tan[so(a2A*) + 3.2]) . (29) 

Again this is invariant under A* — > exp(7r/so)A*. Efimov gave an argument for this functional 
form in Ref. [23]. The formula obtained using Eq. (27) has the same functional form as 
Eq. (29) but does not predict the numbers in Eq. (29). One can also consider the effect of 
corrections to the asymptotic equation of the form 

cos(5oln(p/A,)) Cn .g^s 
cos(soln(po/A*)) {pa2Y ' 

These are subleading in the hmit of large p but become important for p ~ 1/02- The 
coefficient Ci was computed in Ref. [18]. Corrections of this form change the coefficients 
appearing in Eq. (28) but do not change the functional form. 

It would be interesting to know the functional form of r^. An analytic expression for 
the effective range has not been derived. We have extracted from numerical solutions to 
Eq. (1). In Fig. 2, we have plotted r3/a2 as a function of A*a2 for A*a2 between 4.1 and 
92.9. These values of A*a2 differ by 22.7 ~ exp(7r/so), so this interval corresponds to one 
period of the limit cycle, after which rs returns to its original value. In order to compute 
the functional form of rs, it is necessary to know the k dependence of the amplitude K{k, q). 
The asymptotic solution, however, is independent of k. 

Knowing how three-body observables depend on the underlying parameters 02 and A* 
is critical for understanding implications of nonrelativistic scale and conformal symmetry. 
In Ref. [24] , it was conjectured that nonrelativistic systems interacting through short-range 
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FIG. 2. Three-body effective range rs as a function of 02A*. 



forces could be invariant under nonrelativistic scale and conformal transformations in the 
Umit where 02 goes to infinity and the two-body effective range parameters are set to zero. 
This is interesting as it would give an example of a strongly interacting system that is 
nevertheless scale and conformally invariant. The many-body problem in this limit has been 
suggested as a model for neutron matter [25] and could be relevant for gases of trapped 
atoms. In Ref . [24] , it was shown that on- and off-shell two-body Green's functions respect 
the nonrelativistic scale and conformal Ward identities. Ward identities were also derived 
for three-body scattering. Defining the limit 02 — > cxd for the three-body system is subtle 
because of the oscillations in Eq. (22). When we take 02 to infinity continuously, as will 
oscillate rapidly between —00 and -|-oo. While this limit is appropriate for atomic systems 
near Feshbach resonances [7], other hmits may be necessary to find scale and conformally 
invariant theories. One sensible limiting procedure for this purpose would be to take 02 — > 00 
in discrete steps 



lim = lim (02 = expimr / Sr,)a'^) , 



a2— *oo 



(31) 



where n is a sequence of integers and is an arbitrarily chosen starting point for the 
sequence of discrete transformations. An alternative limiting procedure would be to take 
a2 — > 00 continuously while taking A* ^ in such a way that a2A* remains constant. When 
the a2 ^ 00 limit is defined in either of these two ways, observables of the form in Eq. (22) 
live at one point in the limit cycle. 03 does not oscillate when the limit n 00 is taken but 
tends to either zero or ±cxo, depending on the initial choice ag- Obviously obtaining zero 
requires some fine tuning of 02 (or equivalently, A^,). rs will always tend to -|-oo because it is 
positive for all values of A* in its hmit cycle, as is seen in Fig. 2. Higher order terms in the 
effective range expansion will tend to either zero or ibcxD. Because these higher order terms 
can also diverge as 02 —> 00, one needs to understand the behavior of the entire function 
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K[k, k) as 02 oo. It would be interesting to determine this behavior and check if the 
three-body amphtude obeys the scale and conformal Ward identities derived in Ref. [24]. 

One very useful property of the renormalized equations is that it is possible to demon- 
strate analytically that the theory is being renormalized properly rather than having to 
check this numerically with each calculation. The existence of renormalized equations will 
facilitate future numerical work, since there are no delicate cancellations between two dif- 
ferent terms in the kernel. We emphasize that we are not eliminating the three-body force 
but using the renormalization group to simplify the form of the equation and make the de- 
pendence on the three-body force parameter, K^, explicit. In some early treatments of the 
three-body problem [26], the integral equation was solved with a finite cutoff that was tuned 
to fit observed data and then treated as a universal parameter. Though this procedure may 
seem somewhat ad hoc, Eq. (16) shows that this is in fact a rigorous procedure. 

For many applications it will be important to understand how to explicitly renormalize 
higher orders in the EFT. For low-energy neutron-deuteron scattering it is clear that effective 
range corrections are crucial for obtaining accurate agreement with data [1,2,4-6]. Being 
able to renormalize equations analytically at higher orders will facilitate these increasingly 
complex calculations. It would be advantageous to understand the renormalization group 
evolution of higher dimension three-body operators to see if they evolve according to similar 
limit cycles and to understand at what order these enter into calculations. 

A renormalized equation for two-nucleon systems in theories with explicit pions would 
also be of great practical value. Pion exchange gives rise to a potential which must 
be treated nonperturbatively. It is not obvious that the resummed pion graphs can be 
renormalized by a single local operator. If this is possible, one would expect this operator to 
exhibit a limit cycle similar to the one observed in the three-body system, and one would like 
to calculate the evolution of this operator analytically. Investigations along these hues have 
been carried out in a position space treatment using a nonrelativistic Schrodinger equation 
regulated with square well potentials [27]. 

Interesting applications also include coupling photons and weak currents to few micleon 
systems. In these applications one would like to employ gauge invariant regulators such 
as dimensional regularization. This is also important for theories with pions where one 
would like to preserve chiral symmetry explicitly. In this context, it would be useful to have 
renormalized integral equations which are regulated using dimensional regularization. 

We acknowledge useful discussions with P.F. Bedaque, E. Braaten, R.J. Furnstahl, and 
M.B. Wise. This work was supported by the National Science Foundation under Grant No. 
PHY-9800964. 
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